Stochastic behaviour in the edge region of the SINP tokamak 
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Stochasticity is one of the most extensively researched topics in laboratory and space plasmas 
since it has been successful in explaining the various anomalous processes like transport, particle 
heating, particle loss etc. Since there is a growing need for better understanding of this nonlinear 
process, it has led to the development of new and more advanced data analysis techniques. In this 
paper we present an analysis of the floating potential fluctuations which show the existence of a 
stochastic multifractal process along with low dimensional chaos. This has been shown primarily 
by wavelet analysis, and cross checked using other nonlinear techniques. 
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I. INTRODUCTION 

Plasma is a typical complex medium exhibiting a wide 
variety of nonlinear phenomena such as self oscillations, 
chaos, intermittency, etc [H, H, fluctuations in 

the edge region of magnetically confined fusion devices 
have also been associated to nonlinear processes like self- 
organization and chaotic behaviour. Interestingly it has 
also been shown that it is possible to have a coexistence 
of low dimensional chaos and stochastic behaviour 0, [j| . 
The goal of this paper is to present the analysis of float- 
ing potential fluctuations at the edge region of the SINP 
tokamak, wherein these discharges showed an enhanced 
emission of hard x rays signifying the loss of high en- 
ergy runaway electrons [6|. We have deployed several 
techniques, both statistical and spectral to determine the 
nature of stochasticity, and observe the presence of low 
dimensional chaos alongwith stochastic fractal processes 
similar to Bak et al 0, |j| . WTMM technique which has 
been successfully used in other fields, @j H,S EH (3 ^ as 
been used to estimate the multifractal spectrum, and the 
presence of chaos, probably for the first time, in a mag- 
netically confined plasma [lfj. We have cross checked 
these results with other known techniques like nonlinear 
analysis, probability distribution function etc. 

Section |TT] states briefly about the S.I.N. P. tokamak 
and in Section IIIII we have presented the experimental 
results and and discussion and in Section [TVl conclusion. 



II. EXPERIMENTAL SETUP: SINP TOKAMAK 

The experiments were performed in the SINP Tokamak 
(Ro = 30 cm, a = 7.5 cm) which is a small iron core 
machine having a circular cross-section Q . In addition to 
the vertical magnetic field coils it also has an aluminium 
conducting shell 

(i?o = 30 cm, a — 10.9 cm and thickness = 0.7 cm), 
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with four cuts in the toroidal direction and two in the 
poloidal direction respectively. Detail of the SINP toka- 
mak will be found in Ref [6|. The penetration time of 
the conducting shell with cuts (~ 100 [is), keeping con- 
stant the toroidal magnetic field Bt at 0.8 T, toroidal 
electric field (Et) at 30.6 V/m, filling pressure (pfui) 
at 0.2 ± 0.05 mTorr and a at 7.5 cm, and B v was var- 
ied from about 54 mT to 6.75 mT. The edge plasma 
fluctuations are measured using a set of Langmuir and 
magnetic probes and the data were been acquired using 
NICOLET data acquisition system with a sampling rate 
of 1 MHz. In the present work, we report on the analy- 
sis of the floating potential signals from an electrostatic 
Langmuir probe, mounted from the bottom port of the 
toroidal chamber, at r = 7.5 cm. 



III. RESULTS AND DISCUSSION 

An interesting behaviour of the plasma discharges was 
observed in the discharge duration as the equilibrium 
vertical magnetic field, B v was lowered. Fig. [T] shows 
that the discharge duration is almost constant upto 
B v w 20 mT and then increases for B v < 17.6 mT. 
This extension is also clear from the discharge current at 
B v as 13.5 mT [[Figs[2fe)] as compare to current duration 
for 54 mT [Figs EJa)]- The extension in plasma current 
duration was observed after an initial fall to about half its 
peak value. The instant at which the current extension is 
observed to begin is denoted as point B (pt. B) [Fig [2]. 
The horizontal shift in the plasma position (A^ OT .) has 
been shown in Fig. [He) and (g) where, +ve implies an 
outward shift. 

Fig Efb) shows that a few hard X-ray burst were ob- 
served at 54 mT when no extension in the discharge cur- 
rent was observed. On the other hand, the extension 
of the discharge after pt. B is observed to be accom- 
panied by enhanced hard X-ray (HX) bursts [Fig [2{f )] 
which is indicative of loss of highly energetic particles 
from the edge. A characteristic feature of the period of 
extension in these range of discharges is the reduction in 
the electrostatic Langmuir probe floating potential fluc- 
tuations (5<pf) [Fig. 0(h)]. This correlation between the 
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FIG. 3: (color online) Power spectrum, in logarithmic scale 
(a) before pt. B, (b) after pt. B. 



FIG. 1: Discharge current duration (rdisch) vs. B v plot. 
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FIG. 2: (a) and (e) discharge current (I p ) evolution; (b) and 

(f) 3" by 3" Nal(Tl) limiter bremsstrahlung bursts; (c) and 

(g) is the horizontal plasma positioning; (d) and (h) electro- 
static probe floating potential (4>f) signals; for B v = 54 mT 
and B v = 13.5 mT respectively. 



reduction of fluctuations levels with the enhancement of 
bursts of runaway electrons, was a motivation to study 
these electrostatic floating potential fluctuations from a 
time-resolved statistical analysis point of view. 

It is generally accepted that if the power spectrum 
[P(f)) of a signal, obtained from Fast Fourier Transform 
(FFT), decays as P(f) ~ / _ ' 9 where, j3 and / are the the 
spectral index and frequency respectively, then the sig- 
nal shows a stochastic fractal or self-similar behaviour [J]. 
A typical plot of the power spectrum in log-log scale is 
shown in Fig. [31 for the discharge at 13.5 mT and it is 



clear that for the fluctuations before [Fig. [3]Ja)] and after 
pt. B [Fig. 13 a)] follow the power law behavior, which 
indicates the presence of stochastic fractal processes. As 
the power spectrum cannot extract the information re- 
garding the time-frequency simultaneously, the presence 
of sharp transitions and small scale features contained in 
the signal, we introduce more advanced techniques like 
wavelet analysis etc. 
Wavelet analysis: 

Wavelet analysis [H, [3 [l5j], provides a way of an- 
alyzing the local behaviour of functions and correct 
characterization of time series in the presence of non- 
stationarity like global or local trends or biases. One of 
the main aspects of the Wavelet analysis is of great ad- 
vantage is the ability to reveal the hierarchy of (singular) 
features, including the scaling behaviour fig . 

The wavelet transform of a function <f>(t) is then given 
by: 



W*(s,r) = / c(>(t)^ s (t-T)dt 



(1) 



where, <f>{t) is the signal and *£(i) is an oscillating func- 
tions that decays rapidly with time and are termed as 
wavelets, s and r are the scale and time respectively. 
Fig. 2] represents the time-frequency contour plot of the 
power spectrum (| FF^(s,t)| 2 ) obtained from the wavelet 
analysis [Eqn. [Tj, for the discharge at 13.5 mT. For sim- 
plicity the scales in the y-axis have been converted to 
psuedo-frequency (F s ) which has been estimated from 
the relation F s = ^ [l7|, where F c is center frequency 
of the analyzing wavelet and A is the sampling period 
of the signal. The presence of chaos or periodicity can 
be studied using ridge plots obtained from the wavelet 
transform spectrum, which has been discussed in detail 
by Chandre et al [lfj . Fig. 01a) shows the typical ridge 
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FIG. 4: (color onZine^Time-frequency contour plot of power 
spectrum using wavelet transform, intensity of the power have 
given by color bar: (a) before pt. B, (b) after pt. B. Black 
dashed lines show the horizontal ridges connecting maxima 
values in contour plot. 




FIG. 5: Typical multifractal singularity spectrum before 
(solid line) and after pt. B(dot dot line) at 13.5 mT. 



plot at 13.5 mT before pt B. It shows that the most of 
the power is concentrated almost at a constant time scale 
a. The connected horizontal ridges suggest that the elec- 
trostatic floating potential signals are quasi periodic with 
resonance transitions occurring at regular intervals [lol ]. 
For the signal after pt. B, the power is concentrated in 
two or three modes simultaneously at any given instant 
of time 0(b)] indicating the presence of chaos [10j . 

The singularity spectrum, f(a) vs. a, where /(a) is 
the distribution of the singularity strength a, has been 
estimated using Wavelet Transform Modulus Maxima 
(WTMM) method 0, H using the following canoni- 
cal equations, 



«->o in s , 

ln|W*(Mi0))| 
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The singularity spectrum for before and after pt B, for 
discharge at 13.5 mT is shown in Fig. [5] The spectrum 
seems to be slightly asymmetric before pt. B, whereas it 
is almost symmetric after pt. B. The symmetry gives one 
an indication of multiplicative process [HI and hence the 
fluctuations in the extended phase is associated to some 
avalanche phenomena. 

The characteristic of the signals can also be described 
by degree of the multifractality {(3) which is defined by 
the difference between the maximum singularity strength 
(ctmnx) and the minimum singularity strength (a m i n ) [19j, 
M, HE HI- Fig. 1 shows the range of [3 for B v < 15 mT 
is 0.6 — 1.2 for datasets after pt. B and it is 1.2 — 2.7 for 
datasets before pt. B. The decreasing trend in (3 for the 
extended phases indicates that the system has a tendency 
to go towards a stochastic state. 

We estimated the fractal dimension (Dp) and corre- 
lation dimension (D corr = D(q = 2) = 2a — f(a)) from 
the singularity spectrum [5j. Fig.[7Ua) and Fig. [7(b) show 
that in the extended discharge Dp and D corr are in the 
range of 0.86 — 0.84 and 1.8 — 2.0 respectively, indicating 
the presence of complex nature in the signal. 

A crosscheck of the above results for the presence of 
chaos or complexity, can be made by estimating the cor- 
relation dimension (D cr ) and Lyapunov exponent (Al). 
D cr and Xl have been estimated using the Grassberger- 
Procaccia techniques [H, [24[ and the Wolf algorithm 
[1^,[2^| respectively. D cr and Al before and after pt. B 
have been presented in Fig. [5(a) and [5(b) respectively. 
From Fig. [5(a), it is clear that correlation dimension 
obtained using multifractal analysis and Grassberger- 
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FIG. 6: Degree of multifractality (/3) pts before [solid square] 
and after pt. B [open circle] for the various discharges of B v . 
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FIG. 8: (color online)(a) D cr before [black, solid square] 
and after pt. B [orange, open circle] and (b) Al before [black, 
solid square] and after pt. B [orange, open circle] for the 
various discharges of B v . 
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FIG. 7: (color online)(a.) Df before [black, solid square] and 
after pt. B [black, open circle] and (b) D corr before [black, 
solid square] and after pt. B [orange, open circle] for the 
various B v . 



Procaccia techniques are of same order. Fig. [5Jb) shows 
Al is more positive for By < 16 mT after pt. B indi- 
cating chaos. Though we have estimated these exponent 
using insufficient less of data points, the results agrees 
well with wavelet analysis. 

In order to validate our nonlinear analysis we did a sur- 
rogate analysis of the extended discharge regime of Fig. 
[5Ji). The surrogate data has been generated by Phase 
Shuffled surrogate method, in which the phases are ran- 



domized by shuffling the Fourier phases [27|, |28|, l29| , and 
hence the power spectrum (linear structure) is preserved, 
but the nonlinear structures are destroyed [291 ]. D cr 
has been estimated for both the original [Fig[9^b)i] and 
the corresponding surrogate data [Fig E{b)ii], shown in 
FigEJa) by solid circle and open circle respectively. The 
D cr for the original data saturates at higher m, whereas 
in the case of the surrogate data one finds D cr keeps on 
increasing with m. Hence the estimated D cr and Al are 
from nonlinear effects in the system. 

Probabilistic descriptions such as Probability 
Distribution Function [PDF] are at the heart of 
the characterization of turbulence [3l|. Fig. [TOTa) and 
[TDttyshow the PDF at 13.5 mT before and after pt. 
B respectively. The corresponding gaussian fitting is 
shown by dashed line. Both plots show that the PDFs 
are non-Gaussian in nature. Skewness (S) and Kurtosis 
(K) which are measure of nongaussianity are shown in 
Fig. [TTJ f° r the extended discharge which also indicate 
deviation from gaussianity (3lj |. From above analysis it 
is clear that during the extended discharge phase neither 
the system is purely stochastic in nature nor chaotic, 
rather a mixture of both is present. 

The observations of the enhanced energy levels in the 
HX bursts in the extended phase could be a result of loss 
of the high energy particles which are probably generated 
in this phase. The observations of HX bursts can be 
correlated to some growing modes in the dB/dt signals, 
as during the time instants of 1920 /is— 2050 fis, 2090 fis— 
2230 fis and 2320 fis - 2500 fis (Fig. [12]). Subsequently, 
one could infer that some instability could be triggering 
the deconfinement of the particles, which are thereafter 
lost, through some stochastic process at the edge. 
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FIG. 9: (a). Correlation dimension (D cr ) against embed- 
ding dimension for original solid circle and surrogate open 
circle dataset. The dataset used is of the data after pt. B in 
Fig. [2ji) . (b) . (i) Original data and (ii) corresponding phase 
shuffled surrogate data. 
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FIG. 11: Skewness (S) [black open circle] and (K-3) [black, 
solid square] (a) before and (b) after pt. B for various values 
of B v . The dotted vertical line demarcates either side of B v = 
17.6 mT 
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FIG. 10: (color online) PDF of edge electrostatic floating po- 
tentials: before [(a)] and after pt. B [(b))] at B v 13.5 mT. 
Corresponding gaussian fitting is shown by dashed line. 



FIG. 12: Plot 
Plasma Current, 
signals. 



of discharge characteristics depicting (a) 
I p in kA, (b) HX bursts in MeV, (c) dB/dt 



To cross-check whether the discharge could sustain 
any such beam-plasma interactions, we considered the 
conditions that need to be satisfied 13011 : ui re > w 



v e ff > v e , and v beam > 3v c 



where, lu c 



Ve-, Vbeam and v cr are the electron cyclotron frequency, 
electron plasma frequency, electron collision frequency, 
beam velocity and the critical velocity respectively and 



y eff 



'iru> n 



e 3 / 2 A r , A r being the primary run- 



away flux generation factor and e = Eo/E cr , E cr be- 
ing the critical electric field for runaway generation and 
E = Vi oop /2irR Q . Vi oop is the loop voltage. Using 
the experimental results, for beam energy « 1 MeV, 
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n e w 6 x 10 18 m- 3 , B T w 0.8 T, V/ oop w 40 V and 
Z e y^ = 1, we have uo ce rs 140.5 GHz and (x> pe < 100 
GHz and ^ e = 300 kHz and v e fj =1.2 MHz. Hence 
first and second condition are satisfied in the extended 
discharge phase. For the same experimental conditions, 
3v cr [<jJce/^pe] ~ 3 x 10 7 m/s before pt B which cor- 
responds to a beam energy of 2keV and after pt B, 
3v cr [u) ce /ujp e ] « 10 s m/s (corresponding beam energy is 
30 keV). Thus in the current extension phase the more 
energetic electrons will satisfy the third condition. Hence 
it could imply that the loss of runway electrons observed 
in the extended phase could be a result of the participa- 
tion of the higher energy electrons in beam-plasma insta- 
bilities within the plasma column. The edge stochastic 
behaviour possibly leads to the ejection of these runaway 
electrons. 



IV. CONCLUSION 

In the extended discharges of the SINP tokamak, where 
enhanced HX emission were observed, we have shown the 
presence of combination of both of stochasticity and low 



dimensional chaos using various techniques like wavelet 
analysis and other nonlinear techniques like the estima- 
tion of correlation dimension, Lyapunov exponent and 
PDF. One still needs to look into other edge fluctuation 
behaviour, such as the magnetic, density and tempera- 
ture fluctuations, in order to understand the role of the 
stochastic behaviour with the discharge extension, using 
the wavelet transform especially the WTMM method. 
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